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1 Introduction 

Let f(X) e Z[X] be an irreducible polynomial of degree d. It is conjectured 
that, for any integer k > 2, the polynomial f(n) takes infinitely many k-th 
power free values, providing that / satisfies the obviously necessary congru- 
ence conditions. Thus for every prime p we need to assume that there is at 
least one integer n p for which p k \ f(n p ). This problem appears to become 
harder as the degree d increases, but easier as k increases. Thus in 1933 Ricci 
[TO] handled the case k > d, and even proved an asymptotic formula 



N fik {x) ~ A(f, k)x (x^oo) (1) 

where 

Nf,k(x) ■= #{™ e N : n <x, f(n) fc-free}. 
Here the constant C(f, k) is given as 

C(f,k) := - Pf(p h )p- k ) 
v 

where 

Pf(d) := #{n mod d : d \ f(n)}. 

Further progress was made twenty years later by Erdos [2] who showed that 
one could obtain /c-free values for k = d — 1, as soon as d > 3. For such 
k the asymptotic formula ([JJ) was later obtained by Hooley [S]. The next 
development was due to Nair [9] who established jTJ for k > (y/2 - \)d. In 
particular Nair's result shows that k = d — 2 is admissible for d > 24. The 
author [U Theorem 16] then showed how the "determinant method" could 
be applied to the problem, and demonstrated that the asymptotic formula 
remained valid for k > (3d + 2) /4, so that one may take k = d — 2 providing 
only that d > 10. Indeed using methods of Salberger (to appear) one can 
replace these inequalities by k > (3d + 1)/4 and d > 9 respectively. 
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In this paper we show that further progress is possible for irreducible 
polynomials of the form f(X) = X d + c. For these we establish the following 
result. 

Theorem 1 Let f(X) = X d + c G Z[X] be an irreducible polynomial, and 
suppose that k > (5d + 3)/9. Then there is a constant 5(d) such that 

N f;k (x) = C(f, Qx + Oix 1 - 6 ®). 

The implied constant may depend on f and k. 

For comparison with the earlier results we point out that this will allow 
k = d — 2 as soon as d > 6. The result of Erdos handles the case of cubic 
polynomials taking square-free values, and the most interesting open ques- 
tion then concerns quartic polynomials taking square-free values. We would 
therefore like to handle k = d — 2 for d = 4, and one can track our progress 
towards this goal through the historical discussion above. 

There is a related question concerning powerfree values of / at prime 
arguments. Here there is a natural condition that for every prime p there 
should be an integer n p , coprime to p, and such that p k \ f(n p ). With this 
in mind one defines 

N'f k (x) := #{p prime : p < x, f(p) k-iree} 

and 

v 

where 

p'f(d) := #{n mod d : g.c.d.(n, d) = 1, d \ f(n)}. 
The corresponding asymptotic formula 

N' f>k (x)~C'(f,k)ir(x) (x^oo) 

has been proved for k = d by Uchiyama [11] . by a method that also handles 
the case k > d. However it remains an open problem to establish this in the 
case k — d — 1 considered for the previous problem by Erdos and Hooley. 
None the less, important progress has been made by Helfgott [B] and [7], 
showing in particular that the asymptotic formula holds in the case k = 2 
and d = 3. 

Our methods are sufficiently robust that they apply immediately to pow- 
erfree values of f(p). We have the following result. 
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Theorem 2 Let f{X) = X d + c G Z[X] be an irreducible polynomial, and 
suppose that k > (5d + 3)/9. Suppose that for every prime p there is an 
integer n p , coprime to p, and such that p k \ f(n p ) . Then for any fixed A > 
we have 

N' f>k (x) = C'(f,k)n(x) + A (x(logx)- A ). 
In particular this holds for k = d — 1 and every d > 3. 

The preliminary manoeuvres for these problems are straightforward. We 
shall fix the polynomial / (and hence also d) throughout, so that all order 
constants may depend tacitly on / and d. The key fact we shall use is that 



6*|/(n) 



1, f{n) is fc-free, 
0, otherwise. 



It follows that 
with 

and similarly that 
with 



b 

N(b,x) = #{n<x:b k \f(n)}, 
N' f , k (x) = J2Kb)N'(b,x) 



N'(b,x) = #{p<x:b h \f(p)}. 

Clearly N'(b,x) = N(b,x) = for b ^> x d l k . If we denote the solutions to 
f{n) = mod b k by ni, . . . , n r , where r = Pf(b k ), then 

N(b,x) = ^ #{n < x : n = Hj mod b k } 

i<r 

= J>^ + 0(1)) 

i<r 

= xb- k Pf (b k ) + 0( Pf (b k )), 
and similarly, providing that b < (\ogx) 2A we have 

N'(b, x) = E <x:p = Ui mod b k } 



i<r 



E ir(x;b k ,ni) 

i<r, (iii,b)=l 

^p' f (b k ) + A ( Pf (b k )x(logx)-^), 



by the Siegel-Walfisz Theorem. Now, for any £ > we find that 

£ m(6)JV(6, *) = ^%<^> + 0(£ «(*))■ 

&<£ &<§ fc<£ 

The function p/ is multiplicative, with p(p fc ) <C 1, whence 

p(b k ) « e & £ (2) 

for any e > and any square- free b. If k > 2 it follows on taking e — 1/2 
that 

£ = £ + °(£ = <w *) + °<r I/2 ) 

&<§ 6=1 f»£ 

and 

£p# fc )«e 3/2 . 

In particular if we set £ = x 1 ^ 2 we see that 

P(b)N{b, x) = C{f, k)x + 0(x 3/4 ). 

In precisely the same way, if we take £ = (logx) 2A , then 
v(b)N'(b, x) = C'(f, fe)Tr(x) + A (x(log 



xr A ) 



We now consider the range £ < 6 < x 1 71 , where rj is a small positive 
constant. Here we have 

£ rtb)N(b,x)< Yl E (^+°( 1 ))/ ? /( &fc )- 

If we use the bound (jSJ) with £ = |?7 < | this yields 

Kb)N(b, x) « xr 1/2 + x 1 -^ 2 < x 1 -"/ 2 . 

£<&<x 1 -'' 

This bound is satisfactory for Theorem [TJ Since N'(b,x) < N(b,x) we will 
get exactly the same bound in the estimation of N'^ k (x), and again this is 
satisfactory for Theorem [21 
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To complete the proof of the two theorems it will now be enough to show 
that 

x 1 -^<b<^x d / k 

for some 5 > 0, providing that 77 is small enough. By a suitable dyadic 
subdivision we then see that it will suffice to establish the following estimate. 

Lemma 1 Let f(X) = X d + c G 7L\X\ be an irreducible polynomial. For any 
N,A,B EN define 

F{N;A,B) := 

#{(n, a, b) e N 3 : /(n) = ab k , N < n < 2N, A < a < 2A, B < b < 2B}. 

Then if (5d + 3)/9 < k < d — 1 there is a constant S depending on d such 
that 

F(N; A, B) </ N 1 - 5 

forB> N 1 - 5 . 

We have now reduced our problem to one of counting solutions to a Dio- 
phantine equation f{n) = ab k , inside a suitable box. A general procedure for 
such questions is provided by the "determinant method" developed in the 
author's paper [3]. The efficiency of the method depends on the dimension 
of the associated algebraic variety. For f(n) = ab k we are counting integer 
points on an affine surface. Thus far we have made no use of the special 
shape of the polynomial /, but if we observe that f{n) = n d + 0(1) we see 
that (n, a, b) lies close to the weighted projective curve X d = X\X k , where 
Xq and X2 are given weight 1, and X\ has weight d — k. Thus the particular 
form of the polynomial / allows us to consider points close to a curve, rather 
that points on a surface. Reducing the dimension in this way is the key to 
our saving. The procedure is discussed in more detail in the author's work 
[5], to which the interested reader should be directed. 

2 The Determinant Method 

Since f(n) = n d + 0(1) we will have 

N d B -k < A < N d B -k 

for large N . Moreover, since a > 1 we may assume that B k N , and 
indeed we shall assume that 

N 1 ^ < B < N d/k (4) 
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for some positive constant 77. We will choose a parameter K > 1 having 

1 « r ^— « 1, 5 
log N 

and divide the available range for n/b into 0(K) subintervals 

/ = (m N/BK, (mo + 1)N/BK] 

with endpoints defined by integers m in the range 

K < m < K (6) 

We use Fj(N; A, B) to denote the corresponding contribution to F(N; A, B). 
Since /(n) = n d + 0(1) we have n d = ab k + 0(1) and 

(n/b) d = a/b d - k + 0{B- d ). 

It will be convenient to put k = d — j so that 

(n/b) d = a/b> +0{B- d ). (7) 

We now begin the determinant method by listing the points (n r ,a r ,b r ) 
contributing to Fi(N; A, B). Thus the index r runs from 1 to 

R := Fj(N;A,B). 

We choose an integer parameter D > 1 and consider the monomials 

m{n, a, b) = n u a v b w 

for which u + jv + w = D. Thus we may consider D as the weighted degree 
of the monomial, where the variables (n, a, b) are given weights (1, j, 1). The 
number of such monomials will be 

H:= (D-jv + l) = ^r + 0(D) (8) 

v<D/j 

and we label them as mi(n, a, 6), ... , m H (n, a, 6). We now proceed to consider 
the Rx H matrix M say, whose (r,h) entry is mh(n r ,a r ,b r ). The strategy 
of the determinant method is to show that M has rank strictly less than 
H, if the parameters K and D are suitably chosen. If this can be achieved, 
there will be a non-zero integer vector c such that Mc = 0. This vector will 
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depend on the interval /, that is to say it will depend on m . It provides the 
coefficients of a weighted homogeneous polynomial 

Cj(n, a, b) = c h m h (n, a, b) 

h 

such that 

C I (n r ,a r ,b r ) = Q, (r<R). (9) 

If R < H the matrix M automatically has rank less than H. Otherwise 
it suffices to show that any H x H sub-determinant vanishes, and it will 
be enough to consider the determinant formed from the first H rows of M, 
which we shall denote by A. Clearly A is an integer, and our strategy is to 
show that |A| < 1 so that A must vanish. 

We proceed to divide the r-th row of A by b^B~ D for each r < D, and 
similarly to divide the column corresponding to the monomial n u a v b w by 
N U A V B W . Since 

b\ D fnB\ u faBi 



n u a v b w = - — It-t) N u A v B w 
\Bj \bN J \b3AJ 

for u + jv + w = D, this produces a new determinant Ai whose entries are 
of the form rrih(nB /bN, aB^ /WA, 1). Moreover we have 

|A| = | Ai| Y[(b r /B) D Yl N U A V B W < 2 HD P\A 1 \, (10) 

r<H u,v,w 

where 

P = Yl N U A V B W . 

u+jv+w=D 

If we write B = N 13 then we have log A = (d - k/3) log + 0(1), by ©. It 
follows that 

logP = (log N) Yl (u + v(d-kf3)+wf3) + D (l) 

u+jv+w=D 

= (log iV) + (d - kfir 1 + P) + 0(D 2 )} + D (1). (11) 

We now write 

n r B m a r B j N d /m \ d 

MV=^ + S - and ^U=A&\-K +S r) +t - 

Since n r /b r e (m N/BK , (m + 1)N/BK] it follows that 

s r < K~ x . 
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Moreover and (j7|) yield 

a,.!', 1 N d (n r B 



UA AB k \b r N 

and hence 



+ 0(N~ d ), 



Thus the (r, h) entry of Ai will be a polynomial 

/ fc (s r , t r ) = (moK- 1 + s r ) u (N d A- 1 B- k (m K- 1 + s r ) d + t r ) v . 

Clearly f h may depend on h, mo,K,D and d, but it is independent of r. 
Moreover the degree of fh will be at most dD. It follows from (j^J) and (jOJ) 
that N d A~ 1 B~ k <C 1 and moK' 1 <C 1, whence we have the bound | \fh\ \ <Cd 1 
for the height of fh.. 

In order to estimate the size of Ai we will use Lemma 3 of the author's 
work [5J. For each of the monomials s u t v we write 

\\ s u t v \ \ = K- u N~ dv , 

and we list them in order as mi, . . . ,rriH with > > • • •• Then 

according to [5J Lemma 3] we have 

H H 

Ai < D (max ||M|) H JJlKH < c JJlKH- ( 12 ) 

h=l h=l 

To proceed further we shall write K = N K , and note that 1 < k < 1, by 
(jHJ). If we now write m(X), say, for the number of monomials m r = s u t v with 
H^rll > N~ x , then 

A 2 

m(X) = #{(«, !i)6Z 2 : u, u > 0, ku + cfc < A} = — - + 0(A) + 0(1). 

2nd 

If = N~ x ° then m(Ao) > H, while for any e > we will have 

m(Ao-e) <H-1. 
We may therefore deduce that 

A = v / 2«de : +0(l). 

We then find that 

niKii^ - " 



fe=i 



with 



= (™ + dv) + 0(\l) + 0(l) 



KU+dv<\Q 

K +o(xl) + o(i) 



2>nd 
3 

In view of (jHJ), f|10|) . (TXTT) and ffl2l) we may now conclude that 

< ^(l + (rf-A;/3) J - 1 +/3)-^(^) 1 /^3/ 2+OD((logivr i )+0(D 2 ) 
log JM bj 3 

Thus © yields 

l,,dAI < U^ + (d-mr l + ^)- 2 ^dfi\2j)-^ 



D 3 logN ~ Qj x v ^ 3 

+ D ((\ogN)- 1 ) + 0(D- 1 ). 



* = i(l + — + /?) +„. (13) 



We therefore choose 



Ad V j 

with the same small constant rj as in (JH). Then ([5]) will be satisfied, and we 
will have 

log I A I 



< 



D 3 \ogN 

providing that we first choose D = D(f,d,i]) sufficiently large, and then 
ensure that N is sufficiently large in terms of /, d and rj. 

We therefore deduce that A = when K = N K . With this choice the 
matrix M introduced at the beginning of the section will have rank strictly 
less than H, so that all solutions (n r ,a r ,b r ) counted by Fj(N; A, B) satisfy 
the auxiliary equation ©. 



3 Completion of the Proof 

We now complete our estimation of Fj(N; A, B) by considering how many 
triples (n,a,b) can satisfy both the original equation /(n) = ab k and the 
additional equation (j^J). The procedure here will follow precisely that used in 
the author's paper jH §5.3]. Since Cj is homogeneous with exponent weights 
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1) any factor would have to be similarly weighted-homogeneous. It 
follows in particular that Cj(x,y,z) cannot have a factor in common with 
f(x) — yz k . As in [U pages 84 and 85] we find that either 

Fj{N, A, B) < £ (1 + N/B)N £ (14) 

or that there is an irreducible polynomial Gi(X,Y) G Z[X, Y], with degree 
bounded in terms of d and s, but at least d, such that 

(7/(71, 6) = (15) 

for every triple (n, a, b) counted by Fi(N, A, B). 
For a given interval I we will have 

n/bel = (m N/BK, (m + l)N/BK] 

It therefore follows that 



m N 

n T^TT® 

BK 



2N 

<—, B<b<2B. (16) 
K 



It will be convenient to define a linear mapping T : IR 2 — y M 2 by 
rp ( K(2Ny 1 x 1 - {2B)- 1 m Q x 2 

and to consider the lattice 

A = {Tx : x G Z 2 } 

of determinant K^NB)^ 1 . Then if x = (n,b) satisfies (fT6|) we produce a 
point Tx = (ai, a 2 ) G A falling in the square 

S 1 = {(«i,a 2 ) : max(|ai|, \a 2 \) < 1}. 

Let g*- 1 ) be the shortest non-zero vector in the lattice and g*- 2 ^ the shortest 
vector not parallel to g . These vectors will form a basis for A. Moreover we 
have Aig (1) + A 2 g (2) G S only when |Ai| < |g (1) | _1 and |A 2 | < |g (2) | _1 - These 
constraints may be written in the form |Aj| < Lj, for appropriate bounds 
L U L 2 . Since |g( 2 )| > |gW| and |g (1) |.|g (2) | < det(A) < K(NB)- 1 we will 
have L x > L 2 and L X L 2 > NBK~ X . We now write h {i) = T~ x g^ for i = 1, 2. 

These vectors will then be a basis for Z 2 , and if x = Aih 1 - 1 ^ + A2h^ 2 ^ is in the 
region ffTBT) then we will have |Aj| < L, for % = 1, 2. This allows us to make a 
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change of basis, replacing (x 1: x 2 ) by (Ai,A 2 ) so that our constraints on n, b 
are replaced by the conditions |Aj| < L^. 

We therefore proceed to substitute Ai,A2 for n, b in (fl5|) . We may then 
use the bound of Bombieri and Pila [H Theorem 5] to show that the number 
of possible pairs Ai, A2 is <C e max(Li, L 2 ) l / d+e <C e L^ d+e , since the degree of 
Gi is at least d. Thus 

Fj(N,A,B) < £ L\ /d+£ . 

The number L\ depends on the interval /, which is determined by m . We 
therefore write L\ = Li(mo) accordingly. In view of the alternative (fbil) we 
then see that 

F(N, A, B) <C e K(l + N/ B)N £ + ^ £iK)) 1/d+e , (17) 

the range for m being given by ([6]). 

We proceed to investigate the number of choices for tuq which produce 
a value Li(m ) lying in a given dyadic interval (L,2L] say. In the notation 
above, if (n, b) = (xi,x 2 ) corresponds to then 

and L\x 2 <C B. Moreover we will have g.c.d.(xi, x 2 ) = 1. Thus the number 
of intervals I for which L < L\ < 2L is at most the number of triples 
(xi,x 2 ,m ) G Z 3 with g.c.d.(xi,x 2 ) = 1, for which 

x x - ^^x 2 ) < — , Lx 2 < B, and K < m Q < K (18) 

We proceed to consider whether the value x 2 = can occur. If £2 = the 
first of the conditions above would yield Lx\ <C N/K. However we cannot 
have x\ = x 2 = 0, so that we must have L <C N/K whenever x 2 = 0. We 
now recall that L\ 3> L 2 and that L\L 2 ^> NBK^ 1 , whence 

L 2 > NBR- 1 . (19) 

It follows that if x 2 = then (N/K) 2 > L 2 > NBK~ l and hence that 
BK N . However, since .fT = iV K with k given by (fl3l) . we see from (j3J) 
that BK/N tends to infinity with AT, which ensures that the case x 2 = 
cannot arise. 

We now see in particular that the second condition of (ITS1) yields L <C -B. 
If we rewrite the first of the conditions (I18p to say that 

moX2 = N- l BK Xl + 0(BL- 1 ) 
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we then see that each choice for x\ restricts the product m^Xi to an interval 
of length <C B/L, with B/L ^> 1. Moreover m^x^ is never zero. Thus a 
divisor function estimate shows that there are £ (N £ BL~ 1 ) possible pairs 
(x2,m ) for each value of x\. The conditions (1TB]) show that x± <C N/L, so 
that Xi takes 0(1 + N/L) values. This allows us to conclude that the number 
of integers for m which produce a value Li(m ) in the range L < L\ < 2L 
is £ ((1 + N/L)N £ BL- 1 ). 

We can now feed this information into (TTT1) . using a dyadic subdivision 
for the values of Li(rao) to obtain 

F(N, A, B) < £ K(l + N/B)N £ + ^ L 1/d+£ (l + N/L)N £ BL-\ 

L 

in which L runs over powers of 2, subject to the condition L (NBK^ 1 ) 1 ^ 2 
given by (fl9|) . It then follows that 

F(N, A, B) < £ K(l + N/B)N £ + Ll /d+£ (l + N/L )N £ BL \ 

where Lq := max{l , (NBK^ 1 ) 1 ^ 2 }. On taking e = i] we deduce from (j3J) 
that 

F(JV, A, B) N 2r >{K + Ll /d+v (l + N/L )BL 1 }. 
We proceed to analyse our estimate for F(N, A, B) by defining 

if d — kt 



Ad V j 
and q(t) = p(t) + 1 — t. Then 

This is clearly negative if k > j and < t < d/k. Hence if k > d/2 we have 

q (t) > ,Wk) = ^ > o 

for < t < d/k. It therefore follows that KN < B, and hence that Lq < N 
for the relevant range of B. Our estimate now simplifies to give 

F(N, A, B) < r , N 2r '{K + L 2+1/d+r 'NB}. 

This will be of order iV 1-T? if rj > is sufficiently small, and 

sup p(t) < 1 and sup Q(t) < 

l<t<d/k l<t<d/k 

12 



for 

~, N ( A l + t-k(t) 

To handle the condition on p(t) we note that the function attains its supre- 
mum at either t — 1 ot t — d/k. Moreover if v — k/d satisfies 5/9 < v < 1 
we find that p(l) = 9(1 - v)/4 < 1 and 

This latter function is decreasing with respect to v, and takes the value 
196/225 < 1 at v — 5/9. It follows that the supremum is strictly less that 1 
if 5/9 < k/d < 1. 

To verify the condition on Q(t) we note that if 1 < t < d/k then 



dj 2 { 2d\ j J \ j 
1 / 1 \ 2k -d f d-kt 



2d V 2dJ 2d \ j 
1 / 1 \ 2k - d ( d 
~ 2d~\ ~2d) 2d [ + k 



< 
for k > d/2. Thus 

Q(t) < Q(l) 

which is strictly negative for 



S l( 1 -L).( 1 .l 

Ad\ 2d \ d 



. 4d 2d - 2 

J < 



9 2d-l 
This condition is equivalent to 

10d 2 - d _ 5d + 2 2 
> 18d- 9 ~ 9^ + 18d- 9' 

Thus it is necessary and sufficient that 

~ 9 

This completes the proof of Lemma [H and hence also of our two theorems. 
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